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ABSTRACT 

High-energy neutrinos and photons are complementary messengers, probing violent 
astrophysical processes and structural evolution of the Universe. X-ray and neutrino 
observations jointly constrain conditions in active galactic nuclei (AGN) jets: their 
baryonic and leptonic contents, and particle production effi ciency. Testing two stan¬ 
dard neutrino production models for local source Cen A ([Koers fe Tinvakovl I2QQ8I : 
iBecker fc Biermanr] 1200 9b we calculate the high-energy neutrino spectra of single 
AGN sources and derive the flux of high-energy neutrinos expected for the current 
epoch. Assuming that accretion determines both X-rays and particle creation, our 
parametric scaling relations predict neutrino yield in various AGN classes. We de¬ 
rive redshift-dependent num ber densities of each class, from Ch an dra and Swift/BAT 
X-ray luminosity functions ([Silverman et al.l 120081: lAiello et al.l l2009h . We integrate 
the neutrino spectrum expected from the cumulative history of AGN (correcting for 
cosmological and source effects, e.g. jet orientation and beaming). Both emission sce¬ 
narios yield neutrino fluxes well above limits set by IceCube (by ~ 4-10 6 x at 1 PeV, 
depending on the assumed jet models for neutrino production). This implies that: (i) 
Cen A might not be a typical neutrino source as commonly assumed; (ii) both neu¬ 
trino production models overestimate the efficiency; (iii) neutrino luminosity scales 
with accretion power differently among AGN classes and hence does not follow X-ray 
luminosity universally; (iv) some AGN are neutrino-quiet (e.g. below a power threshold 
for neutrino production); (v) neutrino and X-ray emission have different duty cycles 
(e.g. jets alternate between baryonic and leptonic flows); or (vi) some combination of 
the above. 

Key words: black hole physics — neutrinos — acceleration of particles — galaxies: 
active — galaxies: jets — X-rays: galaxies. 


1 INTRODUCTION 

Astronomy has relied heavily on photon-based observations. 
Photons participate in electromagnetic interactions, and 
they inevitably suffer absorption and scattering within the 
emitting sources and in the media along the line of sight. 
Neutrinos are neutral, relativistic particles, but, unlike pho¬ 
tons, only interact weakly with matter. As they are practi¬ 
cally unabsorbed and unscattered over a large distance, even 
propagating through very dense media, they can be used to 
probe the physics of systems at distances as far as the edge 
of the observable universe. Neutrinos are therefore comple¬ 
mentary to photons as astrophysical messenger particles. 
Neutrinos can be generated in violent astrophysical en- 
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vironments. Active galactic nuclei (AGN) and the associated 
jets, together with stellar objects such as pulsars, magne- 
tars, supernovae and q-ray bursters, a re identified as poten¬ 
tial sources of high-energy n eutrinos (iBednarek et al .1120051 : 
IWaxmanl 120071 : iBeckerl 1 20081 ") . AGN are the most luminous 
persistent X-ray sources known. At their cores resides a 
massive black hole (with mass M. ~ 10 6 — 10 9 M©), and 
the accretion of material into their central massive black 
hole powers the AGN activities. The accretion process in 
AGN is often accompanied by a material outflow, which 
manifests as relativistic jets at kpc to Mpc scales. Various 
scenarios for high-energy ne utrino production in AGN jets 
have b een proposed (see e .g. Mannheimlll995l: iMiicke et al.l 


1999: Atovan V Dermerl 120031 : Koers fc Tinvakovl ~ 20081 : 


Becker fc*Biermannll2003 ) . The basic mechanism can be un¬ 


derstood as follows. Charged hadrons, such as protons, are 
first accelerated to very high energies inside the jet. A possi- 
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ble acceleration site is at shocks formed in s ide the jet body 
(e.g. iBlandford V Reesl[l973: lHardeelll979l : iBegelman et al.l 
ll984l ; lBiermann V Strittmatterlll987h . The high -energy pro¬ 
tons accelerated by the jet interact with the ambient par¬ 
ticles (e.g. cosmic microwave background (CMB) photons 
or the baryons in the environments), which generates cas¬ 
cades of lighter children particles and subsequent produc¬ 
tion of charged pions ( 7 r— particles). The decay of these 
charged pion s produces the high-en ergy neutrinos (see e.g. 
IBeckerl 120081 : lArgiielles et al.l l2Q10h . Another possible ac¬ 
celeration site is at the jet base, where accretion inflow 
and relativistic outflow interact. The charged hadrons, pre¬ 
sumably p rotons, are accelera ted in sho cks near the accre¬ 
tion disc (IStecker e t al. Il99ll . see also iNellen et al.l Il993l : 
IStecker &; Salamonlll996h . Through proton-photon (py) in¬ 
teractions with the UV and X-ray photon fields from the 
accretion disc, neutrinos are produced through the decay of 
pions. 

In a theoretical perspective, neutrino production is nat¬ 
urally associated with cosmic rays (CR), as high-energy neu¬ 
trinos are products in the decay chain of energetic parti¬ 
cles produced by interactions between CRs and ambient 
material and photons. In the neutral pion decay follow¬ 
ing the proton-proton (pp) and py interactions, 7 -ray emis¬ 
sion is also produced at comparable energies. The remark¬ 
ably detailed observed CR spectru m extends over eight or¬ 
ders of magn itud e in energy (e.g. lDermer"V Menonl EooiSI: 
IKotera fc Olintol 12011 ). following a power-law with two 


clear breaks, and a suppression of flux tow ards the high¬ 
est energies compatible wi th the GZK effect (|Greise ni ll966l: 
IZatsenin fc Kuz’minl Il96fil ) . The transition from a Galac- 
tic origin to an extragalactic origin is commonly assumed 
to occur around 4 x 10 9 GeV, considering the energet¬ 
ics of known Galactic sources and the non-c orrelation be¬ 
tween local sources an d CR events (see e.g. IBeckerl [2008; 
IKotera fe Olintol 1 20 111 ) . The highest energy CR events are 
therefore tracers of the acceleration processes within the 
sources, however due to cosmic magnetic fields, the CR par¬ 
ticles lose directionality. The neutrinos and 7 -rays produced 
within the sources are however not affected, and whereas 
7 -rays attenuate upon interaction with intergalactic media, 
neutrinos reach us virtually unimpeded. The CR spectrum 
hints that their sources might emit energetic neutrinos too. 
Studying neutrinos and 7 -rays from these sources will enable 
an investigation of the accelerating region within the source 
itself. We take this connection a step further and relate the 
CR emission and its derivatives (neutrinos and 7 -rays) to 
the accretion processes driving the AGN jet. 

Since CR, neutrino and 7 -ray emissions are intrinsically 
linked, the CR and neutrino observations are thus com¬ 
plementary. One may naturally consider that the neutrino 
power scales with the 7 -ray power of the AGN sources. While 
this could be possible for individual sources, the reality is 
more complicated when deriving a scaling relation appli¬ 
cable to the whole AGN population collectively or to an 
AGN subclass population from the 7 -ray observations. For 
instance, by assuming that 10 % of the 7 -ray background 
at the MeV energies is due to non-thermal em ission from 
AGN, such as Seyfert galaxies, ISteckerl (2 005 1) obtained a 
flux <Iv ~ 10 “ 18 GeV -1 cm - 2 s - 2 sr -2 at 100 TeV for the /i- 
neutrinos, comparable the current flux limit of 2.06iQ'g x 
10- 18 (FT/100 TeV ) -2 ‘ 06±0 ' 12 GeV -1 cm - 2 s - 2 sr -2 at‘the 


same energy set by the IceCube exp eriment (lAartsen et al.l 
120151 5 . However, the prescription of ISteckei] (|2005h gives a 
neutrino flux density higher by ^1.5—5 than t he current ob¬ 
serve d limit at PeV energies (see fig. 12 in lAartsen et al.l 
I2015T ). Moreover, it is unclear whether or not 10% of the dif¬ 
fused MeV 7 -rays observed in the sky is non-thermal emis¬ 
sion from the Seyfert AGN and their relation to the neu¬ 
trino generation process. It is also uncertain whether AGN 
in Seyfert galaxies are neutrino source s. Since the IceCub e 
detection of TeV to PeV n e utrinos jA artsen et al.l 120131 : 
IlceCube Collaboration! 120131 : lAartsen et al.l 20151 5 , consis¬ 
tent with an extragalactic origin, several studies have at¬ 
tempted to^pmpointthe source class of these neutrino events 
(see e.g. lHalzenll2014| j. Using the photohadronic interaction 
channel the neutrino flux expected in blazars has been found 
to agree with the IceCube events assuming X-ray and 7 - 
ray emission is produced th rough the 7 r°-decay ([Kraufi et al.l 
I2015l 5 . iDermer et al.l (l2014| j investigates the neutrino output 
by via the photohadronic channel, where the CR protons in¬ 
teract with internal or external radiation fields. They find 
that low-luminosity blazars are poor producers of neutrinos, 
whereas 7 -ray bright flat spectrum radio quasar (FSRQ) 
blazars are promising candidates. In our work we explore 
the neutrino production efficiency in AGN populations by 
focusing on the common engine of the AGN power. X-ray 
and CR emission are both driven by the central accretion 
processes, and the latter will result in the emission of high- 
energy neutrinos. 

In spite of decades of intense observational and theo¬ 
retical studies, in particular in the radio and X-ray wave¬ 
bands, there are still many outstanding questions regard¬ 
ing the dynamical and chemical properties of AGN and 
their jets. We are unsure how much mechanical energy is 
stored in a jet for given observed radio and/or X-ray lu¬ 
minosities (see IWillott et al.l [l 999: M erloni V Heinz! 120071 : 


ICattaneo V Bestl 120091: IShabala V Godfrevl l2013lj and 

put the chemical ingredients in AGN jets (see 
Wardle et al.ll 19981 : iGhisellini V Celottill 200 il : 


Bottcher et al 


120131 5 . In the context of neutrino production, 


we need first to know whether AGN jets are predominantly 
baryonic, leptonic, both baryonic and lepto nic, or electro¬ 
magnetic (i.e. Poynting flux-dominated, see iLovelacelIl976l : 
iLvnden-Belll Il996l : I Nakamura et al.l l2008l j . We also need to 
know if baryons and leptons co-exist in the jet flow and if 
AGN jets have alternating duty cycles of baryonic and lep¬ 
tonic flows, analogous to active and dormant phases in terms 
of X-ray and radio emission. 

AGN are an inhomogeneous class of objects with diverse 
observational properties. For instance, they may be radio- 
loud (RL) or radio-quiet (RQ), and their jets may be weak 
and episodic, or span a large spatial scale and continuously. 
There is no guarantee that neutrinos are produced in the 
same manner across all classes of AGN and that all kinds of 
jets in AGN are equally efficient in neutrino emission. 

In this work, we address the above issues in the neutrino 
output in AGN jets, with an objective to set constraints on 
various scenarios of neutrino production in AGN, using a 
multimessenger approach, which combines the information 
obtained from X-ray observations and neutrino experiments. 
We use the X-ray surv ey observa tions of AGN by Chandra 
([Silverman et al.[2008l j and Swift (lAiello et al.ll2009h and de¬ 
rive the populations of various AGN at different cosmologi- 


© 0000 RAS, MNRAS OOO.fTFTTI 



























































































High-energy neutrino fluxes from AGN 3 


cal epochs. We next apply the neutrino production models 
and determine the neutrino emission from individual AGN. 
From this, we compute the energy spectra and the flux limits 
of the neutrinos generated by different AGN classes and ac¬ 
cumulated throughout the history of the Universe. We then 
compare the flux limits to t he detection limit se t by the Ice- 
Cube neutrino experiment (lAartsen et al.ll2014h . which con¬ 
strains the particle content and physical properties of AGN 
jets, and verifies the neutrino production models proposed 
for different AGN classes. 

This paper is organized as follows: Section [2] discusses 
the argument for AGN as candidate neutrino sources and 
the mechanisms leading to neutrino emission in AGN envi¬ 
ronments. Section [3] describes the two hadron-channel mod¬ 
els for neutrino production in AGN. Section [4| outlines the 
two X-ray surveys that we use, along with a description of 
the X-ray luminosity function (XLF) prescriptions and the 
calculations of the various AGN populations we derive from 
this data. Section [5] presents the high-energy neutrino spec¬ 
tra obtained by our calculations, and Section [ 6 ] concludes 
with the astrophysical implications of our results. 

2 NEUTRINOS FROM AGN JETS 

2.1 AGN as candidate neutrino sources 

Whether a charged particle could attain a certain energy 
depends on the durat ion of i ts con finement within the ac¬ 
celeration region. The iHillasl (|l984| j criterion, which states 

Eq ,max = (fBv , (1) 

gives an estimate for the maximum energy Lg 5ma x of a par¬ 
ticle of charge q accelerated in a region with a characteristic 
size r and a magnetic held B. AGN have emerged as can¬ 
didate neutrino sources, due to their energetic nature and 
ability to accelerate charged particles to energies that facil¬ 
itate the production of high-energy neutrinos. Taking that 
their out hows are hadronic, AGN and their jets are among 
a handful extragalactic source types that ht the require¬ 
ments (see Fig. [I]), and hence are potential 10 20 eV neutrino 
sources. 

AGN emit electromagnetic radiation (photons) over a 
broad waveband from radio, sub-mm, IR, optical, UV to X- 
ray and 7 -ray. They are known as strong X-ray emitters, 
and many AGN are discovered by X-ray observations. In 
practice, compact extragalactic sources observed with a per¬ 
sistent X-ray luminosity above Lx ^ 10 42 erg s -1 could be 
safely assumed as an AGN. 

We would expect variations in the neutrino production 
rates from various AGN classes (Table [J). Often AGN are 
categorized into various subtypes according to their obser¬ 
vational properties at particular wavelengths. A c ommon 
divide is luminosity in radio (e.g. lAntonucdl Il993l j. which 
depends on whether an AGN is bright in radio emission. It 
can be seen releasing its energy in two oppositely directed, 
highly collimated, relativistic jets, perpendicular to the ac¬ 
cretion disc (RL AGN), or with no discernible jet structure 
(RQ AGN). The fraction of RL systems i s about 1CU20 per 
cent of the total AGN pop ulation fe.g. iKellermann et al.l 
Il989l : lUrrv fe Padovanil 1 19951 ) , and the RL fraction is esti¬ 
mated to reach up to 50 per cent for quasars measured in 
X-rays (Idella Ceca et al.|[l994| j. 


Table 1. Classification of radio-loud AGN adapted from 
lUrrv &; Padovanil dl995ll . The three classes of AGN are distin¬ 
guished by inclination of the radio jet to our line of sight. RGs are 
at high viewing angles, and consist of low radio-luminous FR Is 
and higher radio-luminous FR IIs. At lower viewing angles, we 
find RLQs, with SSRQ to FSRQ for decreasing viewing angles. 
At the smallest viewing angles, where the jet is directly in our line 
of sight, we categorize blazars, with lower luminosity BL Lacs and 
higher luminosity FSRQs. The RLQs and blazars are all observed 
with beamed luminosities, and there exists a unification scheme 
across the three types (see the text). 


Type: 

RG 

RLQ 

Blazars 

L x [erg s— 1 ]: 

subclasses: 

10 42 - 10 47 
FR I 

FR II 

10 44 - to 48 - 5 
SSRQ 

FSRQ 

10 44 - to 48 - 5 

BL Lac 

FSRQ 


RL systems are particularly important in the study of 
neutrino output in AGN. Their jets must consist of highly 
energetic, charged particles in order to produce their ob¬ 
served radio luminosities. Here, we present a brief review 
of the various RL AGN subclasses. A common classifica¬ 
tion scheme includes t he orientation and brightness distribu¬ 
tion of their jets (e.g. lBarthellll989l : lUrrv fe Padovanil 1 19951 : 
iTadhuntei] l2008l j . Radio galaxies (RGs) are observed when 
the jet has a viewing angle close to 90°. The active nu¬ 
cleus in these systems is fully obscured or partially obscured. 
RGs are generally separated into two Fanaroff-Riley (FR) 
type s, distinguishable by th e strength of their radio emis¬ 
sion (|Fanaroff fe Rilevlll974| j. FR Type I (FR I) galaxies are 
of lower radio luminosity. They usually show a bright jet 
at the centre. FR Type II (FR II) galaxies are more radio- 
luminous. They have relatively faint central jets, but with 
bright termination shocks at the tip of the jet-blown lobes. 

Radio-loud quasars (RLQs) and blazars are unobscured 
systems in which the jets are aligned along our line of sight or 
close to our line of sight. Their emission is therefore relativis- 
tically beamed. Radio quasars are among the very brightest 
and the most distant objects that we observe. They may 
be separated into lobe emission dominated systems with a 
steep radio spectrum (SSRQ) at higher viewing angles, or 
core emission dominated systems with a flat radio spectrum 
(FSRQ) at smaller viewing angles. Thus, the subclasses of 
radio quasars are also distinguished by the jet orientations. 
If the viewing angle is very small, the jet will be directed 
into our line of sight. This occurs in a blazar. The fraction 
of blazars is no more than 5% of the total AGN population. 
Blazars can be categorized into the high-luminosity FSRQs 
and the low-luminosity BL Lacs. In the framew ork of the 
AGN unification model (|Urrv X Padovanil 11995l b these two 
subclasses are intrinsically considered FR type II/RLQ and 
FR I, respectively, where the jets are aligned in our line of 
sight direction. 

2.2 Neutrino production in AGN jets 

In the hadronic scenario, energetic protons are source par¬ 
ticles for neutrino production. Two interaction channels are 
proposed for neutrino production in AGN environments: 
one with protons interacting with ambient photons (pho¬ 
tons from the accretion disc, synchrotron photons emitted 
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7r —> ii + 

U e~ + z>e + Ufx • (4) 

High-energy CRs are also products of the interactions, as 
escaping neutrons could undergo /3-decays (n —> pe“^ e ), 
leading to emission of neutrinos. 

For neutrinos resulting from pion decay, the ratio of 
neutrino flavours at source is (u e : : v T ) = (1:2:0). Due 

to neutrino oscillations as they propagate through space, we 
expect th e detected ra tio at Earth as (y e : : v T ) = (1 : 

1 : 1) (e.g. lBeckedl2008h . We follow this commonly accepted 
picture in our paper. 

However, deviations due to energy dependences on the 
decay rates and the strength of the source magnetic field can 
lead to energy loss of muons before decay (muon damping). 
In this case, the source ratio is lowered to ( u e \ : v T ) = 

(0 : 1 : 0), as the electron neutrinos from the muon decays 
are of much lower energy than the muon neutrinos produced 
through the more energetic pion decays. This gives the de- 
tect ed flavour distribution as (v e : u u : v T ) = (1 : 1.8 : 1.8) 
(see iKashti &; Waxmanll2005l : |Pakvasall2008h . 


Figure 1. Hillas diagram of the sources which are able to confine 
accelerated protons of maximum energies Fp jma x = 10 16 , 10 18 
and 10 20 eV, with contours of various source candidates, adapted 
fro m fig. 8 inlKotera Sz Olintol (1201 lh . AGN and AGN jets meet 
the lHillad <1984 criterion for energetic protons, and therefore are 
strong candidates for the production of high-energy neutrinos. 
The three points denoted (la), (lb) and (2a) refer to the location 
of Cen A on the Hillas plot with different considerations: (la) 
follows the Hillas criterion (Equation for proton confinement 
at r = 1 32 tv , wher e the confinement is assumed according to 
I Reynoso et all (2011); (lb) at the same location, but taking into 
account energy losses which lowers the maximal energy possible 
in the source, (2a) at a location of r re 3000 r g , w here confinement 
occurs according to iBecker &; Biermannl (l2009l j using the Hillas 
criterion. See Section [3] for details. 


in the jet, CMB photons strayed into the jet); and an¬ 
other with protons interacting with other protons within the 
jet or with protons o f the externa l material entrained into 
the i e t flow (see e.g. Eichler 1979; Mannhei m &; Biermannl 


Il989l : iBegelman et al.l Il990l : Miicke et al. 19991 . and refer¬ 
ences therein). In the py channel, pions are produced via 


p + 7 



( 2 ) 


• P + 7T 


such that two-thirds will follow the 7r°-channel, and the re¬ 
maining third will produce charged pions i t~. The pp interac¬ 
tion also leads to pion production, i.e. pp —> {7r°,7 t + , 7r - }. 
Radiation fields are expected to be strong at the base of 
the jet. In this paper, we consider only models with py in¬ 
teractions and leave those with pp interactions to a future 
study. 

Neutral pions will decay to 7 -rays (ir° —> 77 ), however 
the decay of charged pions produces electrons and neutrinos, 


3 HIGH-ENERGY NEUTRINO PRODUCTION 
MODELS 

In this study w e con sider the model proposed by 
Koers_^_Tm^akgv (120081) and the model proposed by 
Becker &; Biermannl (|2009h . In both models, py interaction 


is the dominant source process, and follows the standard pic¬ 
ture of flavour distribution at observation as outlined above. 
The primary protons are accelerated through shocks in the 
jet, with a power-law energy spectrum. The neutrino flux is 
scaled by CR events detected by the Pierre Auger Observa¬ 
tory (PAO), assuming that the events are of AGN origin. 


3.1 The Koers &; Tinyakov (KT) model 

The model bv lKoers fe Tinvakovl (l2008h studies the relation 
between diffuse and point-source neutrino emission, and uses 
the RG nearest to us, Cen A, as a typical source for neutrino 
production (F ig. [ 2 ]). Cen A, which l i es at a distance of about 
3.4 Mpc (e.g. Ide Yaucouleursl 19791: Soria et all Il996l : llsraell 


U e + lie + Vn, 


( 3 ) 


Il998l : lEvans et al.l 2004 : lHarris et al.ll 2010 h . i^an FR I RG. 
As it is so close to us, it is also a well-used target for neu¬ 
trino studies, as observations can be correlated with its lo¬ 
cation. The mo del thus follows a model of neutrino emission 
from Cen A bv ICuoco fe Hannestadl (120081 . CH), and it is in 
turn based on work on a neutrino produ ction model from 
extragalactic sources of hadronic origin by iMannheim et all 
(1200 11. Following the Hillas criterion, Cen A is a fairly good 
candidate for high-energy particle emission, however when 
accounting for energy losses, it falls within t he boundary of 
maxi mum proton energy E P:inax ~ 10 16 eV (|Revnoso et al.l 
B), and is from these considerations not able to produce 
the highest energy particles. The energy loss calculations 
assume an estimate for the mass of the central black hole 
Mm = 10 8 M©, the location of the confinement and acceler¬ 
ation R ~ 132 r g (where r g = GM % c ~ 2 is the gravitational 
radius of the central black hole), maximum proton energy 
E p , m ax = 2 x 10 7 GeV and formulae for the evolution of the 
bulk Lorentz factor and magnetic field along the jet given in 
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Figure 2. Singl e source spectrum for Ce n A, evaluated at a red- 
shift z = 0 (CH: IKoers fc Tinvakovll2008lf . compare d with the u p¬ 
per fl ux limit for Cen A, determined by IceCube (lAartsen et al.l 


El 


evnoso 


et, all ( 120111 . see Fig. [U point (la) denotes the lo¬ 


cation of Cen A solely based on the Hillas criterion, whereas 
point (lb) shows where it lies if energy losses, with the dom¬ 
inant process being p 7 interactions, are taken into account). 

Protons are confined and accelerated (by shocks) near 
the base of the jet, and interact with ambient X-ray photon 
fields (see Fig. [3]). Following pion production, and subse¬ 
quent decays, neutrinos and neutrons are produced and will 
escape the region. Whereas neutrinos escape unhindered, the 
neutr ons will decay to CR prot ons, which would be observ¬ 
able. iKoers X Tinvakovl (f2008 ) use data from the PAO to 
estimate the ultra high energy cosmic ray (UHECR) flux 
from Cen A, and diffuse UHECR flux, above the threshold 
energy of PAO, UcR,th = 5.7 x 10 19 eV. Due to the com¬ 
mon production path of CR protons and neutrinos, and that 
the emission from Cen A is assumed representative for all 
sources, the UHECR flux is used to scale the neutrino flux. 

There are two breaks in the UHECR proton spectrum, 
which are caused by the change in photo-pion production ef¬ 
ficiency for the protons and neutrons. However, these breaks 
are close in energy, and the model therefore assumes a sin¬ 
gle break energy. This energy is determined through the 
7 -ray break energy, FU .br = 200 MeV for Cen A (see 
IKoers Sz Tinvakovl 20081 . and references therein), such that 
E?cR,br = 3 x 10 8 E 7j br- The all-flavour neutrino flux from 
Cen A, <- A , using the UHECR proton flux <Fp enA as a 
normalization, can then be expressed by 


<£> CenA 

^all 


(£,)« 


& 


min 





(5) 


(see IKoers &; Tinvakovl l2008h . where £ u and £ n are the frac¬ 
tion of proton energy that is converted to the neutrino en¬ 


ergy and the neutron energy respectively, and 


£(FU, EcR,br) 


E v 

77^n-UcR,br 


(6) 


Here the ratio of the average energy of neutrinos and neu¬ 
trons is expressed as rj un = (EC)/(EC) = 0.04, and the frac¬ 
tion of the proton energy transferred to either neutrinos 
or neutrons in the initial interactions is given by % 

0.1/0.5 = 0 . 2 . The UHECR flux from Cen A above the 
thres hold energy UcR.th is <f>^ enA (iU) = 5 x 10 - 21 cm - 2 s -1 
(CH; IKoers fc Tinvakovl 120081 ). The break in the neutrino 
spectrum can then be expressed in terms of the UHECR pro¬ 
ton spectrum and the expression for ratio of average energies 
of neutrinos and neutrons, such that E Ui br = r \ vn EcR,br = 
4 x 10 6 GeV. Due to the neutrino mixing from source to 
observed ratios, the muon neutrino spectrum is one third of 
the full neutrino spectrum, and the flux of muon neutrinos 
from Cen A is 


<J)C enA = a [kt] 



min 



(7) 


with a proton power-law spectrum of index a p = 2.7, and 
that of the neutrino spectrum a u = 1.7. The normalization 
factor, 4 KT 1 contains the scaling through the UHECR flux 
and the combination of energy contribution to the neutrinos 
from the initial particle interactions: 


A KT] = 


(q^p f) ( 


3 

3 x 10“ 


1 GeV -1 cm “s 


E t 


CR,th 


Ei 


is, br 


3>p enA {E th ) 


( 8 ) 


We note that a cut-off of the neutrino spectrum (due to 
it being l imited by the maximum proton energy) is not in¬ 
cluded in IKoers fc Tinvakovl (12008 5, as this will occur at the 
far end of the spectrum. We therefore assume a cut-off at an 
energy of U, max = 5 x 10 19 eV, at which the highest energy 
UHECR event is observed. 

To calculate the diffuse all-sky neutrino flux, 
IKoers V Tinvakovl (|2008h use two source models; one 
with AGN luminosity density of no evolution and another 
with strong evolution (|Bovle V Terlevichlll998h . which lead 
to a diffuse flux ~ 200 — 5000 times larger than that of 
the Cen A flux. This implies a strong detection in either 
case when compared to the current experim ental limit for 
neutr ino detection in IceCube (see Fig. [ 2 j lAartsen et al.l 
l2014h . 

In this paper ho wever, we use the Cen A n eutrino spec¬ 
trum as given in the IKoers fc Tinvakovl (|2008h model, and 
convolve with AGN population densities derived from X- 
ray observations. Because Cen A is a typical source in this 
model, we use the ratio between the neutrino and X-ray lu¬ 
minosities as a reference, and scale the neutrino flux for an 
AGN of a given X-ray luminosity by this ratio: 

/-AGN _ x >CenA / n \ 


using the simplest case, where the scaling factor A = 1, and 


c 


CenA 



( 10 ) 


The X-ray photon intensity of Cen A is l£ enA = Lx/ 47 tD 2 = 
3.47 x 10 -lo erg cm -2 s -1 . Dl is the luminosity distance, 
and for a local source it is the same as the measured proper 
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Figure 3. Schematic illustration of the jet structure (not to scale). On the left is the KT model jet, where shocks at 
the base of the jet accelerate protons that subsequently interact with X-ray photons produced in inverse Compton 
processes in the corona. Neutrons and neutrinos escape the confining regions, however the neutrons suffer beta 
decays before leaving the jet, hence producing a population of CR protons, along with additional neutrinos. On the 
right is the BB model jet. At a few thousand gravitational radii, stable shocks accelerate protons that interact with 
the synchrotron photon field produced by relativistic electrons in the jet magnetic field. Neutrinos escape the jet 
in a collimated beam, whereas protons are continually accelerated along the jet, until they escape the jet as CRs. 
The beam of the CR emission is therefore much larger than that of the neutrinos. Hence, UHECRs may be directly 
observed from AGN with greater viewing angles than sources producing point source neutrinos. 


distance. Thus for Cen A, Dl = 3.4 Mpc = 1.049 x 10 25 cm, 
and a measured (2—J_0 k eV) X-ray luminosity Lx enA = 
4.8 x 10 41 erg s -1 ([Evans et alJ 12004 ). As Cen A is a local 
source, we calculate the single source spectrum at redshift 
z = 0. For an AGN with X-ray luminosity Lx° N , we scale 
the spectrum for a single source with the Cen A brightness 
such that 


t AGN 

a[KT,lin] _ 

^ Lx enA ’ 

giving the AGN fluxes generally, 

$AGN = $ CenA^[KT,lin] 


( 11 ) 


( 12 ) 


This linear scaling then reflects our expectation that a bright 
AGN produces a higher rate of neutrinos than a lower lumi¬ 
nosity counterpart. 


After a comparison of photon optical depths, 
iBecker fc Biermannl (|2009 ) find that the photon field that 
dominates the py interactions is of synchrotron jet origin, 
with optical depth r P1 ~ 1 for a bulk Lorentz factor of the 
jet T ~ 10. The frequency of the photon target field is there¬ 
fore taken to be / = 1 GHz. There are no breaks in the neu¬ 
trino spectrum at these energies, as the break for the radio 
synchrotron photons occurs at much lower energies. 

Furthermore, 29 UHECR events observed by PAO ap¬ 
pear to correlate to AGN in the super-galactic plane. The 
BB model therefore considers the UHECR, and hence neu¬ 
trino, emission to originate in blazars and FR I AGN. The 
neutrino spectrum is then normalized by the proton flux. 

The muon neutrino spectrum traces the proton spec¬ 
trum, and has a cut-off at energies determined by the max¬ 
imum energy of the energetic protons, 


3.2 The Becker & Biermann (BB) model 

W e compare the KT m o del p redictions with a second model 
bv lBecker fc Biermannl ([20091 ). Similar to the KT model, the 
BB model considers neutrinos of hadronic origin in AGN 
jets, and the initial seed protons are accelerated in shocks in 
the jet. However, in this model the peak of neutrino emis¬ 
sion occurs further out in the jet, at the first stable shock, 
a distance of r ~ 3000 r g (see Fig. For an additional 
comparison with the KT model, we have in Fig.[l]added the 
location of Cen A if proton acceleration occurs at a location 
r = 3000 r g (not accounting energy losses). 


$l BB] = 4 bb] j B-“*' exp (—. (13) 

\ t/max / 


The normalization is defined in terms of the redshift 

dependency factor, ^/(/?cr, the ratio of the Lorentz fac¬ 
tors, 7 ^ and ycR of the neutrino and proton (CR) emis¬ 
sion respectively; a measure of the optical depth in the 
source photon field, r P7 ; the Auger threshold energy; and 
flux at energies larger than the minimum of the proton, 
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-E^CR,min — 772 p ^ 7CR ^ GgV, 


^[BB] = 


(«p - !) 

( V" ' 

12 

V^CRy 

x r P7 ( 

m-n \ a P 

4 ) 

1.4 x 10" 

' 10 GeV" 


N(> LcR,min) 


(^PA 0) ap_1 


-1 -2 -1 -1 


(14) 


In the lBecker fc Biermanri (2009 ) model, neutrino emission 
needs to originate in blazar-type AGN to be detected, due 
to the beamed emission. A CR point source may on the 
other hand have a jet misaligned with our line of sight, as 
the emission cone of UHECRs is significantly larger than 
that of the neutrinos (see Fig. EH). The redshift dependency 
(ipv/ip cr) accounts for the difference in flux from neutrinos 
and CRs, based on the redshift evolution of their sources. 
Using radio luminosity functions (LFs) for FR I type blazars 
(BL Lac sources) and FR I RGs, a ratio of the received 
emission of neutrinos and CRs, respectively, is estimated. 
The dependence is expressed as 


r = 




dzdL 


1 dT r dU c 
4nDf dL dz ’ 


(15) 


with the radio LF dT r /d L and the comoving volume el¬ 
ement dVfldz. The factor l/AnD^ takes into account the 
decrease of flux suffered for emission from sources at lumi¬ 
nosity distance LA- To find the ratio between the red shift 
dependency of the emission. iBecker &; Bierm ann (l2009h use 
the fla t spectrum radio source LF fro m iDunlop &T" 3 eacockl 
(1 1990 ). and the FR I radio LF given in lWillott et al.~(|200ll ). 
and estimate a value of ^/(^cr & 0 . 1 . 

We will on the other hand use the source densities, as we 
did for the KT model, based on X-ray observations, and thus 
XLFs for the AGN populations. To scale the n eutrino spec¬ 
tra, w e use the disc-jet symbiosis model (iFalcke fc Biermarml 
Il995l ) that relates the radio power of the jet to the disc lu¬ 
minosity, 


-Kjet 

Ldisc 


= ^d-j 


(16) 


adopt ing disc-jet parameter ftd-j = 0.15 from IBecker et all 
(|2005l ). We use two estimates of the scaling between the ra¬ 
dio and total power in the jet, giving two luminosity scaling 
models. In the first case, we use the synchr otron to total jet 
luminosity relation ([Cavagnolo et al.ll 2010 l ) 


Ljet « 5.8 x 10 


,43 


-Lsynch 

v 10 40 erg s -1 


0.7 


erg s' 


-1 


This gives a disc-jet scaling in terms of the synchrotron lu¬ 
minosity, 

^lynch * 2.0 x lO-^Ldisc) 1 ' 43 erg s ” 1 . (18) 

The s econd scaling model is adopted from IBecker et al.l 
([20051 ). and relates the synchrotron luminosity to the disc 
luminosity by 

^synch * 3.3 x 10 ~ 15 (L^) 1 ' 27 erg s - 1 , (19) 

follow ing IFalcke fe Bier main] (Il995l ) and IFalcke et al.l 
(1995]). The scaling model can be introduced to the neutrino 
energy calculations as the synchrotron luminosity of the 
AGN is proportional to the optical depth for py interactions, 


and the optical depth is also proportio nal to the normaliza¬ 
tion factor for the neutrino spectrum (IBecker fc Biermannl 
120091 ). The optical depth then gives the scaling 


[BBi] _ 


and 


<L 


<t> [ BB 21 = 


Lx 


3.87 x 10 44 erg s ” 1 

_Lx_ 

1.04 x 10 43 erg s _1 


( 20 ) 


( 21 ) 


giving the AGN flux 

^ GN = ^L BB] 4> [BBi] . (22) 


We use a jet Lorentz factor, Tj e t = 10 and a jet half opening 
angle of uj \/2 = FT 4 =0.1 rad. The fraction of the luminosity 
of a knot to the total synchrotron luminosity is e = 0 . 1 , 
and the neutrino production occurs at a distance of 3000 
gravitational radii. 


4 EVOLVING AGN POPULATIONS 


We consider AGN number evolutions from two X-ray sur¬ 
veys, selecting X-ray emitting AGN (thus across t he entire 
popu lation, regard less of radio-loud ness; I Silverman et al.l 
120081 ). and blazars (lAiello et al.ll2009l ). The two give us sets 
of widely different samples of AGN, both varying in luminos¬ 
ity and viewing angle, and where the former contains both 
RL and RQ AGN, and the latter contains only a fraction of 
the RL population. This is accounted for when we derive the 
total number of AGN in the Universe, using the prescription 
of the XLF given in these works. 

Neutrino emission may be direc tly scaled wi th 7 - 
rays originating from pi o n decays (e.g. LHa lzen &; Zasl Il997l : 
lAtovan &; Dermeil l200ll : iNeronov X Ribordvl 120091 ) . How¬ 
ever, only an upper limit can be set on the neutrino emission, 
as some fraction of the emitted 7 -rays would originate from 
the upscattering of e.g. internal synchrotron photons. More¬ 
over, the 7 -ray emission originates in blazar sources, with 
highly beamed luminosities, and are therefore suitable for a 
point source study. In this paper we aim to explore the neu¬ 
trino emission from all jetted AGN classes, hence require a 
survey which is insensitive to orientation effects and obscu¬ 
ration of the jet component. As X-rays trace the accretion 
power of the AGN system, a survey in this waveband fits 
the purpose. 

Furthermore, by choosing X-ray surveys as our AGN 
study, there is no need for complementary observations at 
different wavebands, as X-ray emission implies accretion on 
to a SMBH at the AGN centre. We can therefore be con¬ 
fident that all X-ray luminous objects observed with a X- 
ray luminosity abov e Lx = 10 42 erg s -1 are AGN (see e.g. 


_tv 

iTreister et al l 2004 ). 

I Silverman et al.l ([20081 ) measure the hard (2-8 keV) 
XLF of AGN up to z ~ 5. The sample consists of 682 AGN 
in total, with 31 found at redshifts z > 3. They use the 
Chandra multi-wavelength project to detect high redshift 
luminous AGN (Lx > 10 44 erg s _1 ), an d the Chandra deep 
field to cover the lower luminosity range. lAiello et akl ( 2009 ) 
have used 3 yrs of data from the Swift /BAT survey to select 
a complete sample of X-ray blazars to determine the evolu¬ 
tion of blazars in the 15-55 keV band. The sample consists 
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8 I. B. Jacobsen et al. 


of 26 FSRQs and 12 BL Lac objects in a redshift range of 
0.03 < z < 4.0. 

Both AGN population surveys show that the source 
density evolution of AGN depends on luminosity and epoch. 
The number density can therefore be derived using the XLF, 
assuming that the luminosity distribution of the neutrino 
source s are complete and as inferred by I Silverman et al .1 
(2008|) and lAiello et al.l (l2009h . and can be extrapolated to 
redshifts up to z = 10. 


4.1 The X-ray luminosity function 


The differential XLF of a population is a measure of the 
number of objects per comoving volume and unit luminosity 
interval, as a function of X-ray luminosity and redshift, 

d *{L x ,z) _ d 2 jV(L x ,z) 
d log Lx dV c d log Lx ' 1 1 

The present-day XLF can be expressed as a simple power- 
law (| Aiello et al.H200 9). 


d T(Lx, z = 0) 
d log L x 


= A ln(10) 



(24) 


however, observationally there is a break, and with a high 
enough source count, this break can be seen. A double 
power-law (e.g. lUeda et al.l 12003 ) can fit the observational 
data, with the faint and bright end slopes dictated by Yi 
and T 2 , respectively, for luminosities below and above the 
break luminosity L*, such that 


dfr(L x ,z = 0) = [ /Lx\ Tl [Lx \ T2 

d log L x V L* ) V L * ) 


(25) 


The parameters in the XLFs are determined through 
maximum likeliho od routines, using the MINUIT minimiza- 
tion package (see ISilverman et al.l l2008l : lAiello et al.i r2009. 
for details), and here we use the best-fitting values as given 
in the respective papers (see Table [2]). We also note that 


dT(L x , z = 0) 
d log Lx 


= Aln(10) L x 


dT(L x , £ = 0) 
dLx 


(26) 


The evolution of the XLF depends on the chosen model that 
fits the observations best. The base models are the pure lu¬ 
minosity evolution (PLE) and the pure density evolution 
(PDE), however these are not found to represent the ob¬ 
servational data well. Therefore, modified versions of these 
models are used, either extending the form of the luminos¬ 
ity or densit y evolution, or formulating a combi nation of the 
two (see e.g. lUeda et alJ[2Qo3 : lAird et al.ll2010h . 

The blazar population is found to be best described 
in term s of a modified p ure luminosity evolution (MPLE) 
model ([Aiello et alJ 12009 ) on a double power-law present- 
day XLF (Equation l25l) . where the evolution factor is a 
power law with two free parameters, v\ and V 2 , giving a 
general behaviou r with r espect to redshift, of a form first 
given in lWall et al.l (l2008l b 


dT(L x , z) = dT[L x /e(z),Q] 
d log L x d log L x 


(27) 


e(z) = (1 + z f 1+ * 2Z . (28) 


The fitted parameters are summarized in Table [2] 

We also take a closer look at the best-fitting XLFs of 
the subclasses of the blazars, namely the BL Lacs and FS¬ 
RQs. Whereas the FSRQs are modelled similarly to the full 
blazar sample, the BL Lacs are too few in nu mber, so we 
use the best-fitting simple power-law XLF of lAiello et al.l 
( 2009 ) (Equation [24]), with a simple evolution factor, e(z) = 
(1 + z) Vl (see Table [2]). 

ISilverman et al.H2008l ) determined that the best-fitting 
XLF for their sample is the luminosity dependent density 
evolution (LDDE) model, for which the evolution factor, 
e(z, Lx) is a function of both redshift and luminosity. It 
is convolved with the double power law present-day XLF 
(Equation [25} to determine the population density evolution 
as follows: 


dT(L x , z) dT(L x , z = 0) , 
d log Lx d log L x ’ X 


(29) 


The evolution factor is defined in terms of a luminosity de¬ 
pendent redshift cut-off 2 *, which is further determined by 
a power law of Lx, 


e(z,Lx) 

and 


ui+zr 


je(z*(L x )) 

1 + z 

_ 1 + z*(Lx) _ 


[z < 2 * (L x )] 

\z ^ z*(L x )] 

(30) 


z* (Lx) 



(L x < L c ) 


(L x > L c ) 


(31) 


Another five parameters are therefore introduced when 
evolving the XLF; to determine the redshift cut-off, the 
characteristic luminosity L c , the cut-off redshift z c and the 
strength of the redshift cut-off dependence a. In addition, 
the evolution rates prior to and beyond the redshift cut-off 
z * are given by v\ and V 2 , respectively. 

We assume the XLF at lower redshifts can be extrapo¬ 
lated to describe the high-redshift evolution, and as such we 
span the AGN evolution from redshifts 0 < z < 10. We set 
the upper redshift to z = 10, however note that the oldest 
quasa r is found at a redshift z & 7 thus far dMortlock et al.l 
l201ll ). The lower and upper luminosity bounds on the AGN 
populations are carefully determined, particularly for the 
FSRQ population, as the faint end of the FSRQ XLF col¬ 
lapses towards higher luminosities dependent on the redshift 
bin, seen in Fig. [4] This will be discussed further in the next 
section ( 021 ). 

We carry out our calculations assuming the distribution 
of luminosities obtained from these XLFs is a good repre¬ 
sentation of AGN sources - though we note the possibility of 
missing a low-luminosity AGN contribution in the surveys, 
especi ally at high redshifts. This is mentioned in lAiello et al l 
( 20091) , as BAT is not sensitive to low-luminosity and low- 
redshift sources. The faint end might be under representa¬ 
tive of the real population, as indicated by their results and 
earlier radio-selected surveys of blazars. 

We use the XLFs to calculate the number densities, 
over a range of luminosities and redshifts. This enables us 
to study and compare the neutrino contribution predicted 
from AGN of low and high luminosities, and also from low 
to high redshifts. The comoving volume for a flat, matter 
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High-energy neutrino fluxes from AGN 9 


Table 2. Fitted XLF model parameters. Summary of the fitted parameters for each be st-fitting XLF model, as given in the 
relevant papers: mo del with first letter S refers to model from [Silverman et al.l (12008 ): models with first letter A refers to 
the models given in I Aiello et al.1 (120091) . LDDE: luminosity dependent density evolution; MPLE: modified pure luminosity 
evolution, and PLE refers to the pure luminosity e volution. In brackets a re the source population modelled, such that Type 
I/IIs are described by an LDD E model (model B in lSilverman et al.ll2008h . blazars and FSRQs by MPLE models (best-fitting 
models 7 and 10, respectively: I Aiello et al1l2009l) , and BL Lacs by the PLE (best-fitting model 11). 


Model 

z — 0 Parameters 




evolution 

parameters 



(A, log A) (L 

* , log Lsjs) 

(Ti) 

(T 2 ,T 2 ) 

(vi,vi) 

(u 2 ,u 2 ) 

(Zc) 

(logic) (a) 

SLDDE (Type I/II) 

-6.077° 

44.33 

2.15 

1.10 

4.00 

-1.5 

1.9 

44.6 0.317 

AMPLE (Blazar) 

1.379 x 10“ 7 

1.81 6 

-0.87 

2.73 

3.45 e 

—0.25 e 

- 

- 

AMPLE (FSRQ) 

0.175 x 10“ 7 

2.42 i> 

< —50.0 C 

2.49 

3.67 e 

—0.30 e 

- 

- 

APLE (BL Lac) 

0.830 x 10“ 7 

l.O 6 

- 

2.61 d 

—0.79 e 

- 

- 

- 


Notes. a The value represents the logged normalisation constant, log A, as is given in ISilverman et ahl (I2OO8 ). b The value 


represents the unlogged value of the break lu minosity, L*, wh ere the luminosities are all normalised to L 44 = 10 44 erg s — 1 . c In 
these calculations we used 71 = —50.0 fsee lAiello et al.ll2009t ). d The BL Lac XLF model uses the single power law expression 
( Equation I24D . with index Y 2 . e The blazar, FSRQ and BL Lac XLFs assume an evolution defined by the indices v\ and i) 2 . 


dominated cosmology, is measured ([Peacockl 12007 ) a 
3 (Qz + (fi — 2)[Vl + £lz - l]) 2 




Q 4 (l + z) 3 V 1 + Qz 


(32) 


and we use the cosmological prescriptions given in the rel¬ 
evant papers to maintain consistency of each population 
(Ho = 70 km s - 1 Mpc“\ Ov = 0.7, Q M = 0.3). The X- 
ray luminosities are normalized to L 44 = 10 44 erg s -1 in our 
calculations. 

As an aside, it is worth emphasizing that our calcula¬ 
tions are conservatively based only on the well-understood 
AGN populations. A surprising new radio AGN class ‘FR O’ 
was recognized recently: their radio core s resemble FR I 
cores , lacking extend ed radio emission ([Baldi fc CapettH 
120091 : iGhisellinll 1201 lh , and they may outnumber FR I 
sourc es by ~ 3 - ~ 100 times (ISadler et al.ll20l3 : lBaldi et al.l 
l2015h . Why the FR 0 cores fail to drill their jets farther 
out is unclear - perhaps due to youth, intermittency, in¬ 
terstellar medium obstruction (l ike GPS/CSS sources; e.g. 
lO’Deal fl998l : ISaxton e t al. 20 051) . or intrinsically low T or 
slow SMBH spin ( Baldi et al.ll2015l ). Whatever the reason, 
if FR 0 cores turn out to be as ^-bright as FR I cores, then 
their addition would strengthen our constraints on overall 
AGN neutrino production. The limits also tighten in a sim¬ 
ilar manner if, for instance, our Lx cut-off has underesti¬ 
mated a si gnificant contribution from lower luminosity FR I 
RGs (e.g. (Balmaverde et ahl 120061 : lHardcastle et al. 120091 : 
ICaoetti fe Raiterill2015l ). 


4.2 AGN number density distribution 

We compute the evolutionary tracks and luminosity distri¬ 
butions over several cosmological epochs for all four AGN 
subpopulations, integrating the XLFs with respect to lu¬ 
minosity and redshift, respectively. To obtain estimates for 
the full AGN population, we scale the XLF by a correction 
factor ©cf to obtain the number of all AGN within our 
redshift range. For the RG s ample we a c count for those ob¬ 
scured or too faint following IZinn et ahl (1201 lh , and assume 
that obse rved sources are 10% of the total population. How¬ 
ever, the ISilverman et al.l ( 2008 ) survey collects both RL 
and RQ sources, so we assume the RG population accounts 


for 10% of all X-ray detected AGN ([Urrv &; Padovanill 19951 ). 
The correction factor for RGs is therefore ©cf = 1. In the 
case of blazars we correct for misaligned sources, obtaining 
the correction factor as the ratio of the solid angle of a full 
sphere to the solid angle of the jet projection on to this 
sphere. The viewing angle is uj 0 = Coq/ 2 , so that a jet with 
an opening angle of 2 aq / 2 will be not be in our line of sight 
if the viewing angle is larger than the half opening angle. 
Assuming a modest bulk Lorentz factor T = 10, which re¬ 
lates to the half opening angle by u i / 2 = T _1 % 5? 7, these 
misaligned sources imply that 


©cf = 


4tt 4tt 

Q 2 tt (1 — cos(o;i/ 2 )) 

4n = 4 r 2 . 


2*(w? /2 / 2 ) 


(33) 


This gives a correction factor of 400, which agree s with 
estimates of a few hundred, or 2T 2 (I Aiello et Ml 120091 : 
iGhisellini et al.l l2010l : IVolonteri et al.l l201ll ). We calculate 
the AGN number density evolution over cosmological epochs 
by integrating the XLF into bins of X-ray luminosity to give 
the redshift distribution, such that 


d N(z) 
dz 


©CF 


J 'log L 2 
log L\ 


d *(Lx,z) dV c 
d log Lx dz 


d log L x • 


(34) 


The luminosity dependence of the AGN population found 
in the XLF prescriptions motivates a closer look at the lu¬ 
minosity distribution of AGN in bins of redshift. We assess 
how the dominant luminosity class varies with redshift by 
integrating the XLFs over several cosmological epochs in the 
range of X-ray luminosity adopted in our calculations, which 
gives 


diV(logLx) P 2 dT(L x , z) dVc 

d log L x CF J Z1 d log L x dz Z 


(35) 


We use the appropriate luminosity range, with each 
bin spanning an equal size for a consistent comparison. We 
choose the upper and lower bounds by evaluating the max¬ 
imum luminosity of an AGN, according to the Eddington 
luminosity of a given SMBH mass. For AGN, we assume an 
upper mass of M. ~ 10 9 M©, and we find that the max¬ 
imum luminosity should be about 10 47 erg s -1 . Thus for 
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Table 3. AGN space densities. The space densities calculated for the various AGN populations con¬ 
sidered in a redshift range 0 < z < 10 are shown. The luminosity ranges assumed for each population 
are summarised together with the corresponding unbeamed luminsity. For the blazar population we 
determined this intrinsic luminosity to avoid contamination from other beamed luminous sources at 
high redshift (see the text). 



RG 

Blazars 

FSRQ (FR II) 

BL Lac (FR I) 

L x [erg s 4 ]: 

10 42 -10 47 

io 44 -io 48 - 5 

io 46 -io 48 - 5 

10 44 ' 5 -10 48 - 5 

L x [erg s 1 ] : 

- 

io 40 -io 44 - 5 

10 43 - 8 -10 48 - 5 

o 

1 

O 

o 

Space density [Mpc -3 ]: 

~ 1.36 x 10 -4 

~ 2.26 x 10 -4 

~ 6.84 x 10 -6 

~ 5.79 x 1()- 4 


t hu R L derived AGN population, our range follows that of 
I Silverman et a j] (l2008h . spanning six orders of magnitude. 

In the case of blazars, we need to account for the beam¬ 
in g of these o bj ects, as the quoted X-ray luminosities given 
in lAiello et al. 143) are referring to observed luminosities. 
The jets of observed blazars are beamed in our direction, 
hence the X-ray luminosities we record for the sources are 
greatly enhanced by this phenomenon. As the blazar surveys 
probe the deep past of the Universe, the estimation of the 
population size is based on these luminosities, and we there¬ 
fore make a cut at an intrinsic luminosity of 10 40 erg s -1 . 
We will only take those above this luminosity to be actual 
blazar observations, as we otherwise may confuse some of 
those that fall below this luminosity with e.g. X-ray binaries 
and other luminous objects that could also be observed a t 
these redshifts (e.g. ISwartz et all 120041 : [Feng fe Soriall20llh . 
These sources may also show beamed luminosities compa¬ 
rable to the fainter blazars, but whose intrinsic luminosity 
generally is found at around 10 38 erg s -1 . 

We calcul a te th e intrinsic X-ray luminosity following 
lUrrv fc Shafer! ill 9841) . using the relation between observed 
Lx and emitted luminosity Lx for a relativistic jet, 

Lx = <5^Lx , (36) 

where 5 = [F(l — /3 cos(0))] -1 is the jet Doppler factor, ft is 
the velocity in terms of the speed of light, and the Lorentz 
factor r = [1—/3 2 ] -1//2 , and the viewing angle u 0 = T _1 . The 
exponent g gives the enhancement of the luminosity, and for 
a blazar type in which only one jet is seen, g = 3 + <a, where 
a is the spectral index. This exponent is due to relativistic 
aberration, whereby the emission is beamed forward due to 
the relativistic motions of the jet; contraction of the time 
interval, thus we observe more photons per unit time; and 
the blueshifting of photons, as there are a factor S a more 
photons at the observed frequency than at the emitted fre¬ 
quency. It is found that the observed and intrinsic LFs have 
the same slope for high luminosities, however the observed 
LF will flatten towards lower luminosities because it is sen¬ 
sitive to the l ower cut-off and steepness of the Lorentz factor 
distribution (lListerll2003l ). 

We assume a representative value for the spectral in¬ 
dex of sources in a given AGN population. We use, for the 
subsamples BL Lacs, a = 1.5; for FSRQs, a = 0 .6; and the 
total blazar sample we use a = 1.0 (see fig. 2 I Aiello et al.l 
l2009h . This means that the lowest luminosity bound for the 
full blazar-derived population and the FSRQ-derived pop¬ 
ulation is 10 44 erg s -1 , and for BL Lacs a little higher, 
at 10 44 ' 5 erg s -1 . A further consideration is in order for 
the FSRQ-derived population, because the XLF for these 
sources collapses at lower luminosities, as seen in the FSRQ 


luminosity distribution in Fig. [4] We therefore make a lower 
cut for this population at 10 46 erg s -1 . The assumed X-ray 
luminosity ranges for the AGN populations are summarized 
in Table [3] We find that the cuts we have made do not affect 
the total estimated numbers significantly. 

The large-scale space density (between 0 < z < 10) is 
derived for the various populations, given the comoving vol¬ 
ume contai ned within a redshift of z = 10 is V c % 3.5 x 
10 12 Mpc 3 (IWrightl 1200(3) . These are summarized in Table 
[3j and agree wit h local AGN estima tes from UHECR obser¬ 
vations (see e.g. lTakami et al.ll2012l . and references therein). 
We find that the higher luminosity AGN are preferentially 
found at higher redshifts. The blazar population is domi¬ 
nated by FSRQs at high luminosities (Lx > 10 46 erg s -1 ), 
and peaks here, with the brightest objects most numerous in 
the redshift bin between 4 and 5. This peak moves towards 
lower luminosities in more recent epochs. 

The RG population on the other hand, has a higher 
number density in the redshifts between 1 and 3, and simi¬ 
larly to the blazars, its brightest objects are found to peak in 
this range. In the local Universe the fainter RGs dominate, 
and beyond the peak (i.e. z > 2.5), the source density is de¬ 
clining. The BL Lacs are lower luminosity objects, that are 
declining in numbers with higher luminosity, and the pop¬ 
ulation dominates in the near Universe, in redshifts around 
0 < z < 2. 

Furthermore, at earlier cosmological epochs, we find 
the higher luminosity AGN dominate. There is therefore 
a deficiency of bright AGN in the local Universe. Thus, it 
is suggested that the most luminous blazars formed early 
i n the Universe, f ollowed by a quick decrease in density 
(| Aiello et al.l l2009l ) . Though this implies that black holes 
formed early on in the Universe, and that early Universe 
conditions were favourable to the formation of very lumi¬ 
nous AGN, the observational limitations at higher redshifts 
restricts the study of low-luminosity AGN at the same dis¬ 
tances. 

To check our results from these two AGN population 
estimates, we sum the contributions from either the lumi¬ 
nosity bins or redshift bins, and find that the sums agree, 
giving the space densities as quoted in Table [3] 

In view of neutrino output, we would thus expect a high 
production rate from bright quasars from the epoch of z ~ 4, 
and bright AGN at z ~ 2. If neutrino emission scales with 
the X-ray luminosity, these sources should then produce a 
higher rate of neutrinos than their lower luminosity coun¬ 
terparts. However, the fainter sources are more numerous, 
hence, despite a lower rate of neutrino production, the sheer 
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log [Lx / L44] 



Figure 4. AGN redshift and luminosity distributions. The panels on the left are the redshift distributions of RGs, blazars, 
FSRQs and BL Lacs (top to bottom). The panels on the right give the corresponding luminosity distributions. The overall trend is 
that the brighter AGN peak at earlier times, whereas the low-luminosity AGN are abundant at lower redshifts. Furthermore, the 
bright quasars are seen to dominate at higher redshifts, suggesting the density evolution of the brightest AGN was exceeded by 
fainter AGN at a redshift of z ~ 1. See the text for details. 
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number of these sources would imply a considerable contri¬ 
bution at more local redshifts. 


5 HIGH-ENERGY NEUTRINO SPECTRA 

We make use of the neutrino production models described 
in Section [3] combined with the number distribution of neu¬ 
trino sources from the AGN evolution models in Section 
21 to calculate the expected high-energy spectra, assuming 
that AGN are indeed the dominant high-energy neutrino 
machines (and therefore ignoring other possible sources for 
simplicity). 

We note the distributions of AGN both in luminosity 
and redshift, and find that by applying a luminosity scaling, 
we would expect to find the highest contribution of neutrinos 
from low-luminosity AGN at later times, as these sources are 
the most numerous in this epoch. Due to high-luminosity 
AGN dominating at earlier times (i.e. at redshifts z > 4), 
these should be prominent contributors, as their neutrino 
production rate should be considerably higher than in their 
low-luminosity counterparts. 

The luminosity scaling model is a simple relation be¬ 
tween the neutrino luminosity and the X-ray luminosity of 
the source. This assumes the X-rays originate in the accre¬ 
tion disc, as commonly attributed. Observationally, the ra¬ 
diative and kinetic .jet power corre lates linearly with the disc 
luminosity ([Ghisellini et alJl2014h . The KT neutrino output 
is scaled by a ratio of the AGN X-ray luminosity to that of 
the representative source, Cen A. The other scaling model, 
applied to the BB neutrino spectra, links the disc X-ray 
emission to the jet synchrotron emission, which ties the neu¬ 
trino emission to the disc X-ray luminosity indirectly (see 
Section El. Thus, a more powerful AGN would be brighter 
in X-rays. Similarly, a more powerful AGN will have greater 
potential to accelerate particles in its jet, and contributing 
to a higher rate of interactions. This again leads to an ex¬ 
pected higher rate of neutrino production in these jets. We 
can therefore link the X-ray luminosity of the AGN (i.e. a 
direct measure of the accretion power of the AGN disc) to 
the neutrino luminosity (a consequence of the available en¬ 
ergy in the AGN jet to accelerate and accommodate particle 
interactions of energies related to the AGN power). 

We produce a single source neutrino spectrum, follow¬ 
ing the published models outlined in Section [3] and convolve 
this with the AGN data to obtain emission from the entire 
populations. We scale with luminosity to reflect the influence 
that the source power has on the rate of neutrino produc¬ 
tion. We also correct for cosmological effects. The neutrino 
flux spectrum is required in terms of observable emission 
on Earth, and we calculate our spectra in the source frame. 
Hence, we shift our single source to different cosmological 
epochs, such that the emitted energy at source, Lint, is re¬ 
lated to the energy as we would observe it on Earth, E 0 b s 
through 

L int = £ ob s(l + Z ) • (37) 

We carry out the spectral calculations, and relate the inten¬ 
sity received on Earth, / 0 bs, to the intensity calculated at 
source, Ant, 

/obs = /int(l + «r 4 , (38) 


to obtain the neutrino flux expected to be observed on 
Earth, taking into account cosmological effects such as red- 
shift distortions. For a typical source, the intensity of the 
neutrino emission drops significantly with increasing red¬ 
shifts, and a break in the spectrum will move to lower en¬ 
ergies (as demonstrated in the KT model calculations). We 
obtain the neutrino spectra produced in AGN populations 
from various cosmological epochs by 

El® l model] = Az , (39) 

dz 

and produced by AGN of certain luminosities by 

E 2 $ [model] = ^2 $ AGNd^l0^Lx) A 

d log Lx 

The sum of the binned contributions in each case gives the 
total diffuse emission as we would observe it. We measure the 
spectra against the experimental flux limit set by IceCube , 

L 2 4v [GeV cm -2 s _1 sr -1 ] ^ 1.44 x 1(T 8 , (41) 

in the energy range 3.45 x 10 ^ < LGGeVj < 3.66 x 10 7 , 
determined with 1 yr of data ([Aartsen et al.ll2014h . This is 
an estimate of the minimum neutrino flux required for de¬ 
tection, and therefore gives an upper bound on the neutrino 
flux, as IceCube has detected only a few tens of ev ents so far 
(| Aartsen et al. I [20131 : IlceCube Collaboration! l2013l j. 

5.1 Resultant neutrino spectra 

We present a representation of observable neutrino emissions 
originating in various cosmological epochs, or from a range 
of source luminosities. 

5.1.1 KT model spectra 

The resultant energy spectra expected from the KT model 
prescription is shown in Fig. [5] as the sum of contributions 
binned in source luminosity (solid line) or redshift (dashed 
line). The neutrino emission from RGs far exceeds the Ice- 
Cube limit, with the dominant emission coming from local 
sources, z < 1. At lower neutrino energies [E v ^ 10 6 GeV) 
the greatest contribution coincides with the peak of the AGN 
population, i.e. 1 < z < 2. 

The source luminosities dominating the neutrino emis¬ 
sion are between 43.0 < log Lx < 45.0. Low-luminosity 
sources contribute comparably, due to their abundance in 
the near Universe. The brightest sources are few, and thus 
only contribute significantly at lower redshifts, due to prop¬ 
agation effects. 

The energy spectra of neutrinos originating in the blazar 
and FSRQ populations also exceed the IceCube limit. The 
dominant epoch is 2 < z < 4, which coincides with the peak 
activity of these sources. At the highest neutrino energies 
the local epoch dominates, and the dominant contribution 
at lower neutrino energies extends up to a redshift z < 6. 
The low-energy trend is stronger in the spectra with FSRQ 
origin than that of the full blazar population, and occurs at 
energies E u < 10 6 GeV. 

The luminosities of the blazar sources responsible for 
the majority of the neutrino emission are 45.5 < log Lx < 
46.5. The lower luminosity for the FSRQ population is 
log Lx = 46.0, and the neutrino contribution decreases with 
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Figure 5. Predicted observed neutrino emission from various 
AGN source populations, total diffuse emission from contri¬ 
butions binned in redshift (solid line) and luminosity (dotted 
line). The horizontal solid line is the IceC ube 1-yr (IC-59) neu¬ 
trino detection limit (|Aartsen et alJl2014h . The neutrino emis¬ 
sion expected from the RG population (red lines) gives the high¬ 
est detection, followed by emission from blazar-derived sources 
(green lines) and FSRQ-derived sources (blue lines). The emis¬ 
sion from the BL Lac derived population is the only source 
emission we cannot definitively exclude within this model. The 
dash-dotted line cor responds to the Ice Cube best-fitting diffuse 
neutrino spectrum (|Aartsen et al.1 l2015l b where — 2.06 x 
IQ" 18 [E u /( 10 5 GeV)]- 2 - 46 GeV" 1 cm" 2 s" 1 sr" 1 . 


increasing luminosity, following the same trend as the bright¬ 
est blazars (log Lx > 46.0). 

BL Lacs evolve negatively with redshift, and the neu¬ 
trino flux from these sources follows this trend. The domi¬ 
nant redshift contribution is from the local epoch, with the 
flux decreasing with increasing redshift. The source lumi¬ 
nosity contributions are dominated by the low-luminosity 
sources, and the neutrino emission similarly decreases with 
increasing luminosity. 

The X-ray selected BL Lacs are scarce, and the neutrino 
emission produced in these sources falls below the IceCube 
limit at lower neutrino energies. The only source popula¬ 
tion we cannot definitively exclude within the KT model 
prescription is therefore the BL Lacs. 


5.1.2 BB model spectra 

The resultant energy spectra from the BB model prescrip¬ 
tion is shown in Fig. [6], for the two scaling models used. The 
sum of contributions binned in source luminosity (solid line) 
or redshift (dashed line) are in agreement, and the resultant 
spectra from the two scaling models emphasize the impor¬ 
tance of the luminosity scaling. The neutrino spectra from 
the BB2 model are two orders of magnitude greater than 
those from the BB1 model. Overall, the BB model prescrip¬ 


tion produces lower expected total diffuse neutrino emission 
than that of the KT model. 

The RG population is seen to again greatly exceed 
the IceCube limit. The dominant contribution is from local 
sources, decreasing with increasing redshift, and the bulk 
of the emission originates in bright AGN, with luminosities 
44.0 < log Lx < 46.0. Both scaling models agree on that be¬ 
haviour, however the BB1 model favours the brightest AGN 
over the low-luminosity sources, whereas this trend is weaker 
in the BB2 model. 

The effect of the luminosity scaling is significant for the 
neutrino emissions of blazars. The neutrino flux from blazars 
falls below the IceCube limit using the BB1 scaling, with 
the dominant contribution from AGN with 45.5 < log Lx < 
47.0. The BB2 scaling, however, finds comparable contribu¬ 
tions from lower luminosity sources (44.0 < log Lx < 45.5) 
and brighter sources (46.5 < log Lx < 47.0), and the to¬ 
tal diffuse emission exceeds the IceCube limit by about two 
orders of magnitude. 

The epoch enclosing the peak activity of blazars and 
FSRQs results in the majority of the diffuse neutrino flux 
in these sources, 2 < z < 4. As the FSRQs are already as¬ 
sumed to be the upper end of the blazar luminosities, the 
neutrino emission decreases with increasing FSRQ luminos¬ 
ity. However, the flux from the epoch of 4 < z < 6 is slightly 
higher relative to the blazar population, due to the lower 
luminosities tend towards lower redshifts, and are therefore 
not found in the FSRQ population. 

The neutrino spectra from the BL Lac population re¬ 
flect the source evolution, similarly to the KT model spectra. 
As the population evolves negatively with redshift, the dom¬ 
inant neutrino contribution originates at the lowest source 
luminosity, and in the most local sources. The BL Lacs are 
the only population to be accepted as a possible neutrino 
producing source within the model prescriptions considered 
above. 


6 SUMMARY AND CONCLUSION 

We have calculated the XLF, and derived the total numbers, 
for RL AGN and blazars in different luminosity and redshift 
ranges. These AGN demographics are then convolved with 
the neutrino production model to obtain the muon neutrino 
energy spectra as detected on Earth. It is apparent that 
the neutrino emission received from the epoch of AGN peak 
activity is enhanced, which is a consequence of the X-ray 
luminosity scaling. Additionally, the importance of the lu¬ 
minosity scaling model is highlighted by our results (Fig. 
[ 6 ]). 

We test a number of assumptions in these calculations. 
The AGN source distributions are assumed to derive from 
complete sur veys, and the evolution is correctly repre sented 
by the XLFs ([Silverman et al.ll2008l : lA~iello et al.ll2009h . This 
gives us a fair estimate of the AGN populations at different 
redshifts, and correct proportions over the range of X-ray 
luminosities. 

The modest bulk Lorentz factor, T = 10 is assumed typ¬ 
ical for all AGN. It provides a correction factor for blazars 
accounting for misaligned source s, 0cf = 400, which agree s 
with estimates in literature (e.g. lRachen &; Meszaroslll998f ). 
It also determines the beaming to intrinsic source luminos- 
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Figure 6. Predicted observed neutrino emission from various AGN source populations, total diffuse emission from 
contributions binned in redshift (solid li ne) and lu minos ity (dotted line). The horizontal solid line is the IceCube 
1 -yr (IC-59) neutrino detection limit (lAartsen et al.ll2014l ).The neutrino emission expected from the RG population 
(red lines) gives the highest detection, followed by emission from blazar-derived sources (green lines) and FSRQ- 
derived sources (blue lines). The left panel uses the luminosity scaling BB1, which leads to the acceptance of 
all blazars as candidate neutrino sources. The BB2 model excludes all but the BL Lac p opulation. The dash 
dotted line corresponds to the IceCube best-fit diffuse neutrino spectrum (lAartsen et al.ll2015l ), where = 2.06 x 
IQ" 18 [E u /( 10 5 GeV)]- 2 - 46 GeV" * 1 cm" 2 s" 1 sr" 1 . 


ity relation for blazars, which affects in the luminosity scal¬ 
ing of the re sultant neutrin o spectra. Extreme blazars can 
have T > 20 (lMarscherll2009l ) , but < 30 (|Rachen Meszarosl 
Il998h thus future studies would do well in exploring the ef¬ 
fects of varying Lorentz bulk factors of AGN jets on the 
resultant neutrino emission. 

We assume a correct ion for un de tected sources in the 
RG population following IZinn et ID JmD), assuming the 
detected fraction is 10% of the total population. As we want 
to st udy RGs, we assume th e RL fraction of the survey is 
10% (|Urrv V Pad ovariill 19951 ). thus yielding a correction fac¬ 
tor of ©cf- 

Dim to the scarcity of the BL Lacs population 
lAiello et al.l (120091 ) note that they can only fit a single power- 
law LF. This may reflect an unfair representation of the 
neutrino emission from this sub-population of blazars. Our 
neutrino emission calculations are only as good as the source 
density model, and any conclusion drawn with respect to 
neutrino emission from low or high luminosity AGN, may 
not be valid. It is therefore worth improving the survey data 
to lower luminosities to fill in the lower end of the LF, as 
well as towards higher redshifts. It would also be interesting 
to explore the intrinsic XLF of the blazar populations (e.g. 

I Aiello et al.Loij 120141) . 

The modelling of neutrino production in AGN jets re¬ 
lies on the knowledge of particle interaction channels, and 
the branching ratios of the interaction. The production site 
is assumed to be at the base of the AGN jet as the environ¬ 
ment in the vicinity of the black hole core is sufficiently en¬ 
ergetic for high-energy particle interactions to occur. Models 


of high-energy neutrino output from AGN lobes, the torus 
or intergalactic med ia, due to pp or py interactions are also 
worth studying fe.g. lBecker Biermann|[2009l : lReime 31201 ll . 
and references therein). 


The effect of the luminosity scaling model used in each 
case is shown to be crucial for the estimated neutrino spectra 
(Figs. [5] and [6]). We use a model that will favour neutrino 
emission from bright sources, and suppress emission from 
the abundant lower luminosity counterparts. The KT scaling 
model assumes a linear relation between the source X-ray lu¬ 
minosity and neutrino emission, and is normalized to that of 
Cen A, which is assumed to be a typical neutrino producing 
AGN. The BB1 and BB2 scaling models rely o n the relation 
betw e en the jet synchrot r on and total power (IFalcke et al.l 
Il995l : iBecker et al.l 120051 : ICavagnolo et al. 201 C)h. and the 
jet-disc symbiosis model (IFalcke fe Biermannlll995l ). which 
therefore enables a relation between neutrino emission in the 
jet to the X-ray disc luminosity. 


The resultant diffuse neutrino fluxes predicted using 
these models excee d the observational flux limit set by 
lAartsen et alJ ( 2014* ), implying the following: 

(i) Cen A might not be a typical neutrino source as com¬ 
monly assumed . If Cen A is an exceptio nally efficient neu¬ 
trino machine (|Koers Tinvakovl[20Q 8) the scaling of the 
neutrino yield will lead to an overestimated neutrino flux, 
similar to what our calculations show (Figure [5]). 

(ii) The two neutrino flux models we considered have overes¬ 
timated the neutrino production efficiency. The KT and BB 
models are motivated by an observed correlation between 
UHECRs and local AGN (IPierre Auger Collaboration et al] 
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120071 . 20081 1. The models assume a correlation between CR 
and neutrino emissions due to their common production 
path, e.g. through pj interactions. The KT model is based on 
the observation of UHEC R emission originati n g in t he vicin¬ 
ity of Cen A, however iLemoine &; Waxmanl ( 20091 ) discuss 
the possibility of an accidental correlation between Cen A 
and the observed UHECR events. 

(iii) Neutrino luminosity does not universally scale with the 
accretion power for all AGN subclasses, and hence not with 
their X-ray luminosity. This will require a more complex 
class dependence scaling prescription than the simple uni¬ 
versal scaling that we have used here. The three scaling mod¬ 
els we use are linear (KT) or power laws (BB1, BB2). The 
steepest scaling is given by the BB1 model, and is seen to 
suppress the contribution from the abundant lowest lumi¬ 
nosity sources. As these sources are predominantly found in 
the nearby Universe, a suppression will then enhance the 
neutrino contribution from the epoch coinciding with the 
peak activity of bright AGN, at redshifts 2 < z < 4. 

(iv) Some AGN are not neutrino sources. For instance, there 
could be a power threshold only above which charged parti¬ 
cles could be accelerated efficiently and neutrino production 
could occur. The low-luminosity FR Is may not be suffi¬ 
ciently po werful for the acceleration o f particles to energies 
of 10 20 eV (ILemoine &; Waxmanll2QQ9l ). If FR Is are the par¬ 
ent population of BL Lac sources, then this would also apply 
to t his blazar sub c lass. E nergy loss calculations of the Cen A 
jet ([Reynoso et al.ll2oTlh find that Cen A is unable to obtain 
a proton energy exceeding E Pjmax ~ 10 7 GeV. This is sup¬ 
ported b y indications of low er Lorentz factors in FR Is than 
FR II ([Dermer et al.ll2014h . If FR IIs are the parent popu¬ 
lation of FSRQs, then a neutrino correlation with CRs may 
be weak or negli gi ble, as FR IIs are unfavoured as UHECR 
producers ([Kotera et al.l 12010 ). A highly efficient jet envi¬ 
ronment could lead to the UHECR population decaying be¬ 
fore escaping the confinement, hence only neutrinos would 
be observable. 

(v) Neutrino generation and X-ray emission have different 
duty cycles. Jets may have alternating duty cycle of bary- 
onic and lepton flows, or neutrino production could occur 
only during some fraction of the entire X-ray lifetime of the 
AGN. 

(vi) It is a combination of some of the above. 
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